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[57] ABSTRACT 

A system and method for surface rendering and for the 
display and analysis of arbitrary structures within the inte- 
rior of a solid body is presented. The method uses as a 
starting point a voxel element which is known to be crossed 
by a surface of interest. The surface is rendered in a 
topologically consistent way by connecting encountered 
surface hits at each voxel face. A pair of voxel face hits, 
surrounding a vertex value exceeding the surface value are 
used to form a triangle, spun around a third articulating point 
in the interior of the voxel. Voxel faces sharing the surface 
hits are sequentially examined, and a set of approximation 
triangles is formed, to build a triangularization grid. The 
rendered surfaces are interactively displayed and analyzed 
for global structure properties. 

28 Claims, 20 Drawing Sheets 
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SYSTEM AND METHOD FOR SURFACE 
RENDERING OF INTERNAL STRUCTURES 
WITHIN THE INTERIOR OF A SOLID 
OBJECT 

FIELD OF THE INVENTION 
This invention generally relates to a system and method 
for surface rendering and more particularly to a system and 
method for the display and analysis of arbitrary three 
dimensional surface structures. 

BACKGROUND OF THE INVENTION 

In many practical applications it is required to obtain and 
analyze three-dimensional (3D) arrays of data representing 
one or more physical properties at regular grid positions 
within the interior of solid bodies. Typically, the data is 
obtained by methods such as computer axial tomography 
(CAT), x-ray scanning, nuclear magnetic resonance (NMR) 
imaging systems, or by other non-intrusive mechanisms 
such as ultrasound scanning, positron emission tomography 
(PET), emission computer tomography (ECT) and multimo- 
dality imaging (MMI). Each of these techniques produces a 
planar, grid-like array of values for a continuous sequence of 
planes or slices of the solid object, thus providing a three- 
dimensional array of values. A typical application is in 
medicine where the solid object can be a human body or a 
part of it and the objective is to detect and model certain 
internal structures for diagnostic and surgical planning pur- 
poses. The methods to which the present invention relates 
are, however, equally applicable to other natural or artificial 
bodies and fields of study. 

Examples of the physical properties measured at regular 
3D positions include the coefficient of x-ray absorption in 
the CAT case or, for NMR imaging, the spin-spin or the 
spin-lattice relaxation parameters. In each case, the mea- 
sured physical values reflect variations in the composition, 
density or structural characteristics of the underlying physi- 
cal objects, providing knowledge about their internal struc- 
ture. The 3D data array typically corresponds to a plurality 
of measurements of the given physical property, taken at 
regular intervals along an appropriately selected system of 
3D (X, Y, Z) physical coordinates. Each measurement Vxyz 
of the physical property is thus associated with a specific 
array index position in a cubic lattice within the volume in 
study. Acubically adjacent set of eight such positions defines 
an unit cubic volume of the body under study and is known 
in the art as a voxel. 

It is further known in the prior art to utilize such three- 
dimensional arrays of interior physical values to generate 
visual images of the structures within the body on a display 
system. In order to display two-dimensional (2D) images of 
such 3D interior structures, however, it is necessary to locate 
the position of the surface of such structure within the array 
of physical values. This is accomplished by comparing the 
array values to a threshold value or to a range of threshold 
values which correspond to the physical property values 
normally associated with that surface. 

For example, in the case of the human body, bones, brain 
or any other tissues can be characterized by the known range 
of their density values which is compared to the measured 
and stored 3D array values in order to separate out a surface 
model of their structure. The displayed visual images can 
next be used to assist the medical practitioner in providing 
a diagnosis or in the planning of complex surgical proce- 
dures. 

Once the location of the surface of an internal structure is 
determined, this surface is typically illuminated and shaded 
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so as to give the human eye the impression of a natural shape 
and spatial disposition of the surface when it is displayed. 
Normally, to provide such shading the angular direction of 
a vector orthogonal to the surface at each point on the 
5 surface (the surface normal) is compared to the viewing 
angle of the observer. The intensity of shading can be 
adjusted to be proportional to the difference between these 
angles. 

There are several problems associated with the surface 

10 rendering of structures in the interior of a solid body. First, 
it is necessary to accurately determine which of the voxel 
elements of the 3D array have been crossed by the surface 
of interest and at which physical location. Next, it is impor- 
tant to provide a mechanism for displaying the internal 
structure of a three-dimensional body on the two- 
dimensional computer or TV screen. Another important 
problem is that data acquired from physical data collections 
systems is often noisy and is not capable of being directly 
displayed prior to some kind of a filtering operation. 

2() Several solutions have been proposed in the past to 
address these problems. An approach that is frequently used 
is the one described in U.S. Pat. No. 4,710,876 to Cline et 
al. on the so-called "marching cubes". In this method, the 
surface segment intersecting a voxel is approximated by one 

2S of a limited number of standardized plane polygonal sur- 
faces intersecting the voxel. One particular standardized 
surface tiling is selected from a library of available surfaces 
by a vector which represents the binary differences between 
the threshold value for the given surface of interest and the 

30 eight voxel vertex values. The surface-to-voxel intersection 
coordinates, as well as the vector normal to the surface for 
each such standardized polygonal surface tiling can be 
obtained by table look-up techniques. The final surface is 
assembled as a mosaic using all the standardized triangular 

3J tiles. 

The marching cubes method has found significant accep- 
tance in the industry because it is relatively simple to 
implement and is fast. Several modifications have been 
proposed since 1987. For instance, in U.S. Pat. No. 4,729, 

40 098 granted in March 1988 to Cline et al. a variation of the 
marching cubes methods is proposed, in which non-linear 
interpolation is used to more accurately locate the coordi- 
nates of the standardized polygons. 
Another known method for approximating the surface of 

45 an internal structure is the so-called "dividing cubes" 
method, disclosed in U.S. Pat. No. 4,719,585 to Cline et al. 
In this method, the values at the vertices of the voxel are 
used to interpolate in three dimensions values at regularly 
positioned intra-voxel sub-grid locations. The interpolated 

50 sub-grid values can then be used to locate the surface 
position and calculate the normal vector more accurately. 

In order to distinguish between different internal struc- 
tures having the same or similar physical property and thus 
having close or similar values on the three-dimensional grid, 

55 U.S. Pat. No. 4,751,643 to Lorensen discloses a technique 
for labeling surfaces with similar property values and using 
adjacency criteria in respect to a "seed location" in the 
particular structure of interest to segregate the desired sur- 
face from all of the labeled surfaces. 

60 U.S. Pat. No. 4,791,567 to Cline et al. discloses another 
technique of segregating similar internal structures by deter- 
mining connectivity from adjacency information. U.S. Pat. 
No. 4,821,213 to Cline et al. discloses yet another technique 
for differentiating internal structures io which a linear path 

65 is made through the data array to locate and label all 
different structures along the scan line by counting structure 
interfaces. 
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In the methods described above, voxels within the solid 
object under study are typically traversed sequentially and 
progressively, so that every voxel within the data space is 
checked for a crossing by the surface of interest, and is 
subsequently assigned to belong either to the outside or to 
the inside of the surface. Because of the large volumes of 
data to be processed, the methods are computationally 
intensive and relatively slow. An improvement in this 
respect is described in U.S. Pat. No. 5,166,876 to Cline, 
issued in November 1992. The disclosed method signifi- 
cantly speeds up the rendering of a surface within a bound- 
ary by only processing voxels on the surface. This approach 
saves significant computation time but still uses the same 
case table and the limited set of approximating surfaces 
which were described in the original marching cubes patent. 

As noted above, the marching cubes surface rendering 
method is widely accepted but has certain deficiencies. 
These deficiencies are generally related to the fact that the 
fifteen canonical cases used in the surface segments' 
approximation do not cover the entire spectrum of possible 
tiling combinations. This results in occasional approxima- 
tion errors, occurring at lower levels of resolution, specifi- 
cally if a voxel face is intersected by the surface twice 
resulting in four crossings ("hits") with the voxel edge. The 
errors manifest themselves in incorrect surface rendering, 
visually appearing as unexpected spots, tears and missing 
tile elements. The problems are aggravated if the collection 
data is noisy and/or with low contrast. 

A number of corrections to the marching cube algorithm 
have been offered in the past. Notable is a paper by Neilson 
and Hamann, "The Asymptotic Decider: Resolving the 
Ambiguity in Marching Cubes," in the Proceedings of 
Visualization 1991, pp. 83-91. This and other publications 
show that the case table approach produces inconsistent 
results in some cases. Thus, no topologically correct algo- 
rithm for isosurface tiling has been offered so far. 
Accordingly, reliable measurements of internal structure 
parameters, such as surface area and volume, are also not 
available. 

The display of 3D graphic images has been principally 
driven by the goals and directions of computer aided design 
(CAD) and computer aided manufacturing (CAM) systems. 
Systems accepting data in several basic formats have been 
developed for displaying solid bodies and manipulating 
images in various ways to create models of parts being 
manufactured, or internal organs being monitored. 

The present systems for surface rendering and display are, 
however, still incapable of meeting the requirements for a 
fast and truly interactive display of continuously monitored 
3D data. The complexity of the task can be better appreci- 
ated by noting that a typical surface rendering operation 
involves processing on the order of 512x512x512 pixels. In 
the presence of noise the task becomes even more compli- 
cated because of the necessity for additional noise filtering, 
smoothing and data interpolation, which takes processing 
time and often results in the loss of important structural 
details. 

It is thus desirable to provide a method for rendering 
topologically correct surface approximations of arbitrary 
internal structures within a volume of interest and a system 
for the display and analysis of such surface structures. It is 
also desirable to provide high quality, interactive display of 
the surface structures of interest from a set of collected data 
and minimize the effects of possible noise interference. 
Additionally, it is desirable to provide an efficient imple- 
mentation of the method which facilitates the use of spe- 
cialized hardware for real time data processing. 
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SUMMARY OF THE INVENTION 

Accordingly, it is an object of the present invention to 
provide a system and method for topologically correct 
surface rendering of internal structures within a volume of 
5 interest. 

It is another object of the present invention to provide a 
system and method for the display of three-dimensional 
information. 

jo It is also an object of the present invention to provide a 
system and method which can be implemented in special- 
ized electronic hardware. 

It is another object of the present invention to provide a 
graphics imaging system capable of being interactively 

IS used, and providing high quality images with textural shad- 
ing and other visual clues to the user. 

It is yet another object of the present invention to provide 
medical practitioners with the abilities to emulate surgical 
procedures graphically prior to undertaking invasive patient 

20 treatments. 

Additionally, it is an object of the present invention to 
provide a plurality of three-dimensional surface views from 
a single set of collected data. 

A further object of the present invention is to provide a 
25 system for analysis of internal structures within a 3D solid 
body. 

Lastly, but not exclusively it is an object of the present 
invention to minimize the effects of possible noise interfer- 

30 6nC6 - 

In accordance with a preferred embodiment of the present 
invention, a novel method for rendering three-dimensional 
surface structures comprises the following steps. Signal 
patterns representing the value of at least one physical 

35 property associated with regularly spaced grid locations in a 
three dimensional body are stored in a memory means. The 
regularly spaced grid locations define voxel elements within 
the body, such that each voxel is associated with eight vertex 
values and six voxel faces. The signal patterns that represent 

40 values lying within a specified range associated with the 
surface structure of interest are next selected for further 
examination. 

An initial grid location which corresponds to a position 
within or on the surface structure is first selected and the 
45 eight vertex values associated with a voxel at the initial grid 
location are selectively retrieved. A threshold value is sup- 
plied to determine which voxel vertex values are inside and 
which voxel vertex values are outside of the surface struc- 
ture. 

50 For each voxel face having at least one voxel vertex inside 
and at least one voxel vertex outside of the surface structure, 
there must be at least one pair of intersects between the 
edges of the voxel face and the surface structure. The 
coordinate values of these intersects between the voxel face 

55 and the surface structure are then computed. Preferably, the 
computation of these intersect coordinate values is done by 
a linear interpolation of the voxel vertex values. Other 
methods of computation can similarly be used. 
A triangular approximation to the surface structure is then 

60 formed using the pair of intersects between the edges of the 
voxel face and the surface structure as two vertices of the 
triangle with a third vertex in the interior of the voxel. The 
step of forming triangles is repeated for each pair of voxel 
face intersects on the six faces of the voxel. The physical 

65 location of the third vertex of each triangle is identical and 
is computed as the centroid of all interpolated intersect 
coordinate values for the voxel. 
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The method further comprises the steps of identifying all 
adjacent voxels sharing a face intersect with the surface 
structure. Next, for each identified adjacent voxel, the steps 
of retrieving the eight vertex values, determining the inside 
and outside vertex values, computing the coordinates of the 5 
voxel face intersects and forming the triangles is repeatedly 
applied until all face-adjacent voxels within the three- 
dimensional body having a face intersect with the surface 
structure are processed. In the step of computing the inter- 
sect coordinates, only values that have not previously been 
computed are processed. 

The proposed method of surface rendering further com- 
prises the step of displaying a two-dimensional image of the 
approximated surface structure. Preferably, the surface 
structure is displayed by utilizing only the voxels identified 
by the processor as being crossed by the surface structure. 1S 

In a second preferred embodiment of the present 
invention, a method for generating of an internal surface 
structure within a three-dimensional body comprises the 
following steps. Three-dimensional signal patterns repre- 
senting the value of at least one physical property associated 20 
with regularly spaced grid locations in a three dimensional 
body are stored in memory means. The regularly spaced grid 
locations define voxel elements within the body, such that 
each voxel is associated with vertex locations and voxel 
faces, each two intersecting voxel faces defining a voxel 25 
edge. Each voxel edge passes through two voxel vertices. 

The signal patterns which represent values lying within a 
specified range associated with the surface structure of 
interest are next selected. Two voxel vertices of a voxel edge 
for one of said stored voxels are then retrieved; and it is 
determined which voxel vertex of the selected voxel edge is 
inside and which voxel vertex of the selected edge is outside 
of the surface structure. 

For each selected voxel edge having at least one voxel 3J 
vertex inside and at least one voxel vertex outside of the 
surface structure, a coordinate value of a first hit location 
where the surface structure crosses the voxel edge is next 
computed. Preferably, the computation of the intersect coor- 
dinate values is done by a linear interpolation of the voxel w 
vertex values but other computation methods can be used as 
well. 

Next, each of the two intersecting voxel faces that define 
the voxel edge where a hit was identified is examined for a 
second hit location which defines an intersect between the 45 
surface structure and the voxel face. 

For each intersect of a voxel face with the surface of 
interest, a triangle with a vertex in the interior of the voxel 
is formed to obtain an approximation to the surface struc- 
ture. The location of the third vertex of each approximating 50 
triangle is computed as the centroid of the voxel face 
intersects such that the voxel face intersects form a closed 
contour around the inside voxel vertex. 

The method further comprises the steps of identifying all 
adjacent voxels that share a face intersect with the surface 55 
structure. Next, for each identified adjacent voxel, the steps 
of selectively retrieving two voxel vertices of a voxel edge, 
determining the inside and outside vertex values, computing 
coordinates of the first hit locations, examining the voxel 
faces defining the voxel edges, and forming the approxima- 60 
tion triangles are repeatedly applied until all face-adjacent 
voxels within the three-dimensional body having a face 
intersect with the surface structure are processed. In the step 
of computing the intersect coordinates, only values that have 
not previously been computed are processed. $5 

In accordance with the present invention, the method 
further comprises the step of displaying a two-dimensional 
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image of the approximated surface structure. Preferably, the 
surface structure is displayed by utilizing only the voxels 
identified by the processor as being crossed by the surface 
structure. 

In accordance with a preferred embodiment of the present 
invention, the proposed method for 3D surface rendering is 
implemented on a system for analysis and display of three- 
dimensional surface structures which comprises the follow- 
ing elements. First, storage means are provided for storing 
3D signal patterns representing the value of at least one 
physical property associated with a 3D body at regularly 
spaced grid locations within the body. The signal patterns 
are provided by a three-dimensional raw data source, such as 
MRI, CAT, x-ray, ultrasound or others. Aset of 3D-adjacent 
grid locations define a volume element referred to as voxel 
having associated with it voxel edges and voxel faces. 

The system also comprises means for defining a plane of 
pixel positions having preselected values and corresponding 
to one plane of the spaced grid locations. Also provided is 
means for determining in the defined pixel plane, a pixel 
associated with a voxel that is intersected by a surface 
structure of interest within said body. 

The system further comprises processing means for com- 
puting coordinate values of the intersection between the 
surface structure of interest and the intersected voxel edges 
and for computing the coordinate values of at least one 
articulating location positioned in the interior of the inter- 
sected voxel. Second storage means for storing the com- 
puted coordinate values along with connectivity information 
about joining the computed coordinate values into a set of 
triangles is also provided. These triangles approximate the 
intersection between the surface structure of interest and the 
voxel. Each triangle has one edge on a voxel face and one 
vertex at an articulating location in the interior of the 
intersected voxel. Preferably, the processor means further 
comprises generating means for generating normal vectors 
associated with the computed coordinate values. The com- 
puted normal vectors are similarly stored in the second 
storage means. 

In a preferred embodiment of the present invention, the 
processing means of the system further comprises means for 
identifying which of the voxels sharing a voxel face with the 
intersected voxel are intersected by the surface structure of 
interest. To provide for a better data management of the 
system, the processing means further comprises means for 
generating, for each voxel that is intersected, a plurality of 
index values representing information about the shared 
voxel face intersects. These index values are also stored in 
the second storage means and are used to indicate to the 
processor which voxels and voxel faces have been pro- 
cessed. 

The processing means of the system comprises means for 
simultaneously comparing all vertex values of an intersected 
voxel to a threshold value. This comparison is done inde- 
pendently for each voxel face thus enabling a parallelizable, 
high speed implementation of the proposed system. 

In order to facilitate the display of the 3D surface structure 
of interest, the system of the present invention further 
comprises display processor means which receives data 
from the second storage means and converts it into a display 
format. The display format data is supplied to means for 
interactively displaying surface structures of interest. 

In accordance with another embodiment of the present 
invention, the planar arrangement of the regularly spaced 
grid-locations is other than a rectangular array. For example, 
the planar arrangement of the grid locations may have a 
hexagonal, or other pattern suitable for processing shape. 
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The system and method for surface rendering and display 
of the present invention have a number of practical appli- 
cations. For example, the internal data structure representa- 
tion allows the user to simultaneously process and display 
two or more isosurfaces. As a result, the system provides a 
mechanism to count 3D disjoint structures within the interior 
of a body under study by enumerating the independent 
structures being processed. 

The system and method of the present invention can in 
addition be used to compute and display crosssections of 
disjoint surface structures. Such surface structures need not 
only be internal structures of the object under study. While 
the description of the proposed method is presented in terms 
of the measured values of a given physical property of the 
object, other types of surfaces can similarly be computed 
and displayed. For instance, this can include isogradient and 
isocurvature surfaces giving information about the locations 
within the body having identical rate of change in the density 
values or curvature values, respectively. The cross-section of 
such surfaces with the structure surface can be used to 
automatically compute cardinal landmarks on the object 
surface. 

In a different application, the system of the present 
invention can be used to compute global structure param- 
eters such as the total surface area, and the structure volume. 
Information necessary to compute these quantities, includ- 
ing the coordinates of a set of physical points which belong 
to the surface is stored and available for each processed 
structure. The resulting parameter estimates can be used to 
define characteristics of the object under study, such as 
volume occupancy, average structure size and others, which 
characteristics have a wide range of applications, especially 
in 3D pattern classification. 

Lastly, but not exclusively, the system of the present 
invention can be used to observe time changes in the interior 
structures of a given solid object by knitting together "snap- 
shots" of the surface of the structure and displaying these 
snap-shots in a sequence. An exemplary application of such 
system use is the monitoring of a patient's heart beat pattern. 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other advantages, features and objects of the 
present invention are described in the following description 
of the preferred embodiment in which: 

FIG. 1 is a detail of surface tiling errors in a prior art 
surface rendering method. 

FIG. 2 is a block-diagram of a system for carrying out the 
method of the present invention. 

FIG. 3 is an isocontour plot of an MRI axial slice of a 
brain used for plane selection preview. 

FIG. 4 is a three dimensional voxel with elements defined. 

FIGS. 5(A) and (B) are perspective views of a voxel 
neighborhood and voxel adjacencies for three dimensional 
voxel elements, respectively. 

FIG. 6 shows intersection of a voxel by surface structures. 

FIG. 7 illustrates voxel face hits assumptions. 

FIG. 8 shows ambiguity in voxel face contour lines. 

FIG. 9 is a flow-chart of the SpiderWeb A method. 

FIG. 10 illustrates spinning approximation triangles using 
an articulation point. 

FIG. 11 illustrates spinning approximation triangles for a 
voxel face contour ambiguity using SpiderWeb A. 

FIG. 12 is a flow-chart of the SpiderWeb B method. 

FIG. 13 shows a two-dimensional adaptive sampling grid. 
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FIG. 14 shows triangle vertices drawn in a T-mesh order. 
FIG. IS is an approximating surface computed with full 
precision. 

FIG. 16 shows the approximating surface of FIG. 16 using 
5 midway hit location speed processing. 

FIG. 17 shows triangle face and vertex normal vectors. 

FIGS. 18 and 19 illustrate the medical imaging surfaces 
rendered in accordance with the method of the present 
1Q invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

Referring now to FIG. 1 a surface detail is shown which 
illustrates tiling errors in a prior art surface rendering 
method. The errors are observed as bright or very dark spots 
on the surface, indicating a tear in the represented surface or 
a missing approximation tile. The present invention is 
designed to remove such errors by presenting a topologically 
consistent method of building a triangularization grid which 
represents an arbitrary three-dimensional structure. 

Illustratively, the invention is practiced on the raw data 
obtained from non-intrusive measurements of at least one 
physical property within a planar slice through a body under 

2J study. A typical source of the raw data is an MRI or CT 
system, an x-ray, ultrasound or another type of a scanner, 
capable of converting measurements of some physical prop- 
erty of the body, at regular grid-like locations within the 
planar slice, into an array of discrete values. Other non- 

30 intrusive data collection mechanisms such as positron emis- 
sion tomography (PET), emission computer tomography 
(ECT) and multimodality imaging (MMI) systems can also 
be used as a data source. 
The operation of the data source may be carried out for a 

3S sequence of slices of the solid object, thus providing a 
three-dimensional array of values representative of the spa- 
tial structure of the object. Alternatively, it is possible to 
collect data over various periods of time, for the purpose of 
measuring temporal changes in the body under study. 

40 The raw data is stored in a first storage means, preferably 
indexed by a slice number. The storage means preferably has 
a fast access time and may comprise magnetic storage, 
random access memory or other known in the art storage 
media. Advantageously, the storage means also comprises 

45 means for converting the raw data into an application 
specific internal data structure, in order to take advantage of 
the specific hardware configuration used. For example, the 
internal data structure provides the capability of individually 
addressing and retrieving individual pixel elements as well 

50 as sets of spatially adjacent pixel configurations, referred to 
as voxels. 

In accordance with a preferred embodiment of the present 
invention, the first storage means also provides the capabil- 
ity of indexing each vertex, each voxel edge and voxel face, 

55 in order to facilitate the computation process. 

The indexing of the voxels will be better understood from 
FIGS. 4-6. FIG. 4 shows a voxel element consisting of a 
cubic arrangement of eight data pixels and presents element 
definitions consistently used in the following description. 

60 Each pixel element is associated with a triplet of indexes 
(i,j,k) used to determine the location of the pixel in a 
three-dimensional indexing system. By convention, the 
lower left vertex of the front face is selected as a reference 
point. The indexes of the adjacent pixels to the left, to the top 

65 and towards the back plane are shown to increase by the 
length s, which corresponds to the number of display pixels 
inscribed within two consecutive data points. 
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It is to be understood that FIG. 4 only serves as an 
illustration and other possible conventions are possible, as 
long as they consistently applied. Thus, a different vertex of 
a voxel element may be selected as a reference point, and the 
direction of index increase may be selected opposite to the 5 
one shown in FIG. 4. Furthermore, the length s may corre- 
spond to substantially different physical distances depending 
on the particular application. In addition, even though in 
FIG. 4 the length is identically labeled in all three index 
directions, in practice this may not necessarily hold true. The to 
reason is that the data sampling step in each coordinate 
direction depends on the corresponding spatial resolution of 
the raw data source in the specific index direction. 

In accordance with the present invention, voxel elements, 
such as the vertices, edges and voxel faces, illustrated in 15 
FIG. 4, can be referenced by a data processor in groups or 
individually. Label information is associated with each 
voxel element, so that vertex values connected along edges 
of the voxel, or voxel faces are marked as being connected. 
Similarly, other data structures which are used in a group 20 
during the processing are marked as being connected. A 
second storage means associated with the data processor 
comprises a memory which is allocated for each data 
structure, so that the processor is capable of determining 
whether or not a particular element of the data structure has 25 
been previously processed. For instance, each time a specific 
voxel edge is retrieved, a mark is left indicating for all 
associated data structures that this element has been pro- 
cessed. Thus, when the next element of a connected data 
structure is visited, the processor need not revisit the original 30 
voxel edge. 

FIG. 5(A) is a perspective view of a voxel neighborhood. 
FIG. 5(B) depicts possible combinations for voxel pairs, 
which define a vertex adjacency, an edge adjacency and a 
face adjacency, depending on whether only a vertex, an edge 35 
or a voxel face, respectively, are shared. 

In reference to FIG. 4, a voxel element is shown inter- 
sected by two separate isosurfaces. Voxel vertices, having 
values above an isosurface threshold are termed "inside" 4Q 
vertices, while voxel vertices having values less than the 
isosurface threshold are termed "outside" vertices. The 
intersection of an isosurface of interest with a voxel edge is 
referred to as a "hit". Six such hits, indicated as H0-H5 are 
shown in FIG. 6. The intersection between the surface of 4J 
interest and a voxel face is referred to as a "face intersect" 
and, as seen in FIG. 6, connects two voxel edge hits. 

As previously noted, the first storage means comprises 
means for converting the raw data into a specific internal 
data structure, which allows all voxel elements of the data 50 
array to be consistently indexed by the data processor for 
identification purposes. For example, the face intersect 
H0-H1 on the left face of the voxel in FIG. 6 is shared with 
the adjacent left voxel and, in fact, occupies the same 
physical location. By using a consistent indexing scheme, as 55 
well known in the computer graphics art, the data processor 
need only compute the hit locations once, and use their 
values when processing all corresponding edge adjacent 
voxels, as denned in FIG. 5(6). 

In accordance with the present invention, the data pro- 60 
cessor implements a topologically consistent method of 
connecting the isosurface hits on voxel edges into a con- 
tinuous triangularization grid representing the surface of 
interest for analysis and display. The method requires a 
starting point, or a "seed" pixel element on a data plane, 65 
selected by a plane selecting means. Using only local 
information about the adjacent pixels, the method "spins" a 
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web of triangles, closely approximating the surface of inter- 
est. The processing sequence and visual appearance of the 
rendered surface resemble the construction of spider web, 
giving the process its "Spiderweb" name. 

Some assumptions about the topology of the approxi- 
mated surfaces used in the surface rendering method of the 
present invention are illustrated in reference to FIGS. 7 and 
8. In FIG. 7 depicts a voxel face with vertices V0-V3, edges 
EO-E3 and a clockwise scanning direction, as indicated by 
the arrow. A hit HI is shown on voxel edge El. In accor- 
dance with the present invention, it is first assumed that the 
crossing isosurface of interest is sufficiently smooth. This 
assumption imposes a limitation on the surface curvature. 

Secondly, an isosurface is assumed to have no 
singularities, breaks, or holes, so that given a hit on one 
edge, there will be at least one more hit on the same voxel 
face. As a result of the these assumptions, in accordance with 
the present invention hits on a voxel face must come in pairs. 
Referring back to FIG. 7, this implies that there must be at 
least one additional hit on the voxel face. If no such hit, or 
an odd number of hits is found, the processor will indicate 
an error. 

FIG. 8 illustrates another possible situation, when all four 
edges of a voxel face have isosurface hits. As shown, this is 
only the case when diagonally opposite voxel vertices are 
both "inside," or both "outside" the surface of interest, and 
the two diagonal pairs are not on the same side of the 
surface. Lacking global information about the isosurfaces, 
the exact direction of surface intersects cannot be deter- 
mined. This ambiguity results in four possible cases of face 
intersects, indicated as cases A, B, C and D in FIG. 8. Cases 
B and D will be reconsidered in connection with the Spi- 
derWeb A, and SpiderWeb B processing methods described 
next. 

SpiderWeb A 

FIG. 9 is a flow-chart of a first embodiment of the surface 
rendering method in accordance with the present invention, 
referred to as SpiderWeb A In this method, no decision 
about the exact direction of the face intersects in the case of 
four hits is made. This corresponds to case D in FIG. 8, 
where the four hits are connected, thus minimizing the 
probability is that an error is made. The algorithmic steps in 
this case are as follows. 

In step 200 of the method, the pixel elements correspond- 
ing to one voxel are simultaneously retrieved from the data 
storage. The vertex elements of the voxel are next compared 
in step 210 to a preselected threshold which corresponds to 
the surface of interest, and each vertex is assigned as either 
inside or outside with respect to the surface. 

If an intersection of the surface with the voxel is not 
detected, the method proceeds to retrieve a face adjacent 
voxel in accordance with predefined voxel sequencing rules. 
If a surface intersection is detected at step 220, the locations 
of all edge hits of the voxel, which have not been previously 
processed, are computed. Preferably, the position of a hit 
along an edge is computed by using a linear interpolation of 
the values at the two vertices defining the edge. For example, 
if a voxel edge with a hit has an outside vertex at position 
V^j, with a value of V 0 , and an inside vertex at position V,] 
with a value of V 1( and the threshold value is t, so that 
V 0 <t<Vj, the location of the edge hit Ex can be calculated 
using the formula: 

Ex-V M HKy-V^-VMV>-VJ 

Naturally, other methods of estimating the position of the 
hit, such as non-linear interpolation, may be used, as is well 
known in the art. 
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The next step 240 involves the calculation of an Articu- first voxel face and then around the edges of the second 

lating Point (AP) which is located entirely in the interior of voxel face. Hit locations are calculated by linear interpola- 

the voxel. The articulating point in the SpiderWeb A method tion as discussed above. Since face intersects on a voxel face 

is computed as the mean (centroid) of the system of all occur in pairs, the scanning of each face stops when the first 

connected hits in a voxel. 5 additional hit is located. The processor then forms a triangle 

The articulating point is used in the next step 250 as the having one vertex at the intersection point on the edge that 

third vertex of a set of surface approximation triangles, each was originally hit, a second vertex at the location of the first 

triangle being defined by the articulating point and each of additional hit and a third vertex at an articulating point with 

the possible different pairs of the connected hits in the voxel. an initially undetermined physical location inside the voxel. 

As can be appreciated by those skilled in the art, if only to if me foregoing procedure identifies any single hits on 

two edge hits are found for each voxel face, as for example another face of the voxel, that face is likewise scanned to 

shown in FIG. 6, the articulating point will lie on the plane i ocate the second hit; and the two hits on that face and an 

defined by the hits; and the approximation triangles will be articulating point at an undefined location within the voxel 

co-planar. An illustration of the SpiderWeb A two-face-hits constitute another surface approximation triangle. If a pair of 

approximation surface is shown in FIG. 10 where the 15 hits are identified on another face of the voxel, these hits and 

articulating point is located near the center of the voxel at the an articulating point likewise constitute another surface 

centroid of the four hits and four approximation triangles are approximation triangle. This process continues until all 

denned, each with the articulating point as one vertex and surfaces of the voxel are accounted for. The exact location 

with two connected hits on one of the four faces as the other 0 f the articulating point for all triangles of the voxel is then 

two vertices. 20 computed as the centroid of all edge hits, and information 

If, on the other hand, four hits are encountered in a voxel a b ou t the triangles is stored. Triangle and vertex normals, to 

face, the resulting approximation surface triangles will form be discussed next, are computed at step 370. 

a pyramid, touching at its tip (the AP) another approximation Fof each yoxel e(Jge with a surface ^ me pw(xss tepea(s 

surface. An example is shown m FIG. 11. itself in the face adjacent voxels. In particular, in step 380 a 

In step 260 of the method, triangle and vertex normals for 25 se( of ed from (he face adjacent VQXel b examined until 

each of the approximate triangles are computed The &u Me of the isosurface of j^est are 

purpose of these normals will be discussed in more detail processed . FinallV) in step 390 the processor exits to storage 

sol°70 ofXTSm tfFIG P V Pr ° C£S ' means 60 and the dis P lav Pressor 70. 

'"Jthe following Lp 270, for each intersected voxel face, 30 Unu , ke the S P id " Web A method, this embodiment of the 

, - ,. , ... , . j, present invention does not necessarily involve simultaneous 

the face adjacent voxel is retrieved and processed. In accor- r . , . . . , ' ,, . ... , , 

, ... .. . . .. ..... c ii i processing of entire voxels. Instead, after the initial seed 

dance with the present invention, the indices of all voxel • , , . , ■ , .i 

. ... , , , , ij.u. Pixel is specified, the proposed method tracks along the 

faces that have already been processed are marked so that f . ... ,, ..... 

. , „ , , , j j • isosurface of interest in accordance with the adopted index- 
processed voxel faces need not be considered again. . , ..... , , , c r ... ... 

r- ii j . c r. ni . ■ c , . ing scheme, only visitmg voxel edges and faces with hits. 

Finally, data from the SpiderWeb processor is forwarded 35 6 ' 3 6 e . 

to storage and display at step 270. Both SpiderWeb A and B processing methods reflect the 

SpiderWeb A processing method has the disadvantage of uncertainty of the voxel face contour line ambiguity, illus- 

creating "tunnel" voxels and "bubbles" in the case of a voxel m FIG. 8. In the embodiment of the present invention, 

face ambiguity. However, in the case of a noisy data, this usin g SpiderWeb A, no assumption is made about the exact 

property may actually reduce the visual impact of the noise, 40 direction of the face intersects, resulting in a three dimen- 

since its effect is localized to a particular voxel element and sional "bubble". In the second embodiment of the present 

is not visible to the user. invention, using SpiderWeb B, a systematic choice to "sur- 
round inside vertex" is made based on physical consider- 

SpiderWeb B ations. Yet another possibility to avoid the face intersect 

FIG. 12 is a flow-chart of the second embodiment of the 45 ambiguity is to increase the resolution of the space sampling 

surface rendering method in accordance with the present grid until no voxel face has four hits per face. Thus, in 

invention, referred to as SpiderWeb B. In this method, a accordance with another embodiment of the present inven- 

consistent decision is made about the exact direction of the tion the resolution of the sampling grid is adaptively refined, 

face intersects in the case of four hits on a voxel face. The FIG. 13 shows a two-dimensional adaptive sampling in 

decision corresponds to case B in FIG. 8, where the four so which the upper left corner is illustratively subdivided until 

edge hits are connected in a way so as to surround the inside adequate resolution is achieved. 

voxel vertex. The algorithmic steps in this case are as The rectangular array of sampling elements considered 

follows. thus far is only illustrative of the different sampling grids 

In step 300 of the method, a starting or seed vertex is which can be used in accordance with the present invention, 

selected. In the next step 310, the three edges associated with 55 Well known in the literature is for instance the two- 

the seed vertex are scanned so as to retrieve the vertex dimensional hexagonal sampling grid that can reduce the 

elements at the opposite ends of those edges, and in the number of required samples in a plane by about 13.4% (see, 

following step 320 those vertex elements are compared to for example Dudgeon et al., "Multidimensional Digital 

the pre-selected threshold. The hit locations are calculated Signal Processing," 1984, pg. 39-47). Other sampling grids 

by means of a linear interpolation, as discussed above, or 60 can also be used, as known in the art. 

estimated in any other suitable manner, as known in the art. It is important to note that while the processor computes 

Triangles are formed in SpiderWeb B in the following a large number of triangles, information for which needs to 

four steps which are performed for each edge that is hit. be stored, the actual computation process can be signifi- 

First, in step 340, the two voxel faces that define the edge cantly sped up due to the fact that separate voxel faces can 

that has been hit are scanned for additional hits. This is done 65 be processed independently. Thus, the proposed method can 

beginning with the inside vertex of the edge that has been hit be implemented in a parallel architecture resulting in a 

and proceeding from vertex to vertex around the edges of the significant improvement in terms of processing speed. 
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As shown in FIG. 2, illustrative apparatus for the practice 
of the invention comprises a first storage means 20, a plane 
selection means 30, control unit 50, a Spiderweb processor 
100, a second storage means 60, a display processor 70 and 
a structure analysis means 80. 5 

Plane selection means 30 operates under the control of an 
operator from control unit 50 to select data from individual 
data planes stored in data storage 20. The operator may 
preview data planes corresponding to a particular physical 
location of interest in the body under study. Since structures 10 
having different physical properties exhibit different ranges 
of measured values, the operator can isolate a specific 
structure of interest by selecting a particular threshold or a 
range of threshold values. For example, the operator could 
isolate bones from the surrounding soft tissue in medical 1S 
imaging. 

FIG. 3 illustrates an isocontour plot of an MR1 axial slice 
of the brain used for a plane selection preview in accordance 
with a preferred embodiment of the present invention. In this 
embodiment, using selection means 32 in FIG. 2, the opera- 2 o 
tor can choose a specific plane of interest and display it on 
a monitor. A preselect threshold value is supplied to display 
processor 70 and is used to display the locations (the contour 
lines) at which the selected data plane is crossed by a surface 
structure having the same physical property value. The 2 s 
operator can then identify a particular contour line and 
specify the coordinates of a desired starting point ("seed"). 
The coordinates of the seed are next supplied to processor 
100 to initiate the rendering of the desired surface structure. 
If no contour line of the surface of interest is identified, the 30 
operator can interactively specify a different data plane or 
modify the threshold value until a desirable contour is found 
and a seed element selected. 

In an alternative embodiment of the present invention, the 
selection of a starting point to render an internal surface of 35 
interest can be done using computer selection means 34. In 
this embodiment, a threshold value supplied by the operator 
is used in sequentially scanning the data array for elements 
that cross the threshold value. The first pixel having the same 
value as the threshold value corresponding to the surface of 40 
interest is selected as a starting point and its memory address 
is passed to processor 100 which generates an internal 
surface of interest. After completing the given surface, as 
indicated by absence of face-adjacent voxels that share a 
face intersect with voxels being intersected by the surface, 45 
processor 100 sends a signal back to selection means 34 
prompting for selection of the next starting point. Since all 
voxels visited by the processor 100 are "flagged" to indicate 
that they have been processed, selection means 34 only 
identifies the next unprocessed starting point on a surface of so 
interest. Its address is passed to the processor for rendering 
the next surface structure. The process continues until no 
voxel within the stored data array having a vertex value 
corresponding to the desired threshold is left without being 
flagged as processed. This provides an indication that all 55 
disjoined surfaces that satisfy the threshold conditions have 
been processed. 

Storage means 60 is used by the SpiderWeb processor to 
store data about the interconnection structure of the pro- 
cessed isosurface. In accordance with the processing 60 
method, this may include information about the visited 
voxel vertices, edges and faces, as well as the computed 
locations of edge hits and articulating points, and informa- 
tion about the interconnection of different elements into 
individually addressable data structures, such as the approxi- 65 
mation triangles spun from each face intersect. Also stored 
in storage means 60 are the so-called shared normal vectors 
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associated with the articulating points and voxel edge hit 
locations. The shared normal vectors are computed in pro- 
cessor 100 as the average of the set of unit normal vectors 
for each triangle having a vertex at the respective hit or 
articulating point. 

The display of the generated approximation surfaces is an 
important element of the proposed invention. While both 
SpiderWeb A and B methods in general exhibit insensitivity 
to noisy data which is results in a better display, they both 
generate a large number of triangles, thereby increasing the 
computational burden of the system. In a number of cases, 
as discussed earlier, such triangles are in fact co-planar and 
in principle need not be processed separately. The specific 
data structure used in the methods may, however, be used to 
significantly reduce the computational burden at the display 
stage. 

Reference is now made to FIG. 13 in which a "T-mesh" 
ordering of approximation surface triangles is shown. The 
T-mesh run illustrated in the figure is a contiguous strip of 
triangles, ordered in a "zig-zag" fashion. Aset of three initial 
triangle vertices, labeled 0, 1, 2 is first drawn. After the last 
triangle vertex 2 is specified, the triangle 0 is drawn, using 
the connecting edges between vertices as delimiters. Since 
the edge (1-2) is shared by triangles 0 and 1, when rendering 
the latter triangle, one need only specify one more point 
(vertex 3). Thus, after vertex 3 is specified, triangle 1 can be 
drawn. In the next step, vertex 4 builds triangle 2, by only 
specifying one additional point. The process continues, until 
all triangles of the T-mesh run are drawn. At each step only 
one additional triangle vertex is specified and the display 
processor never needs to draw an edge twice. 

The SpiderWeb method of surface rendering results in a 
data structure especially suitable for display using the T-run 
approach illustrated in FIG. 13. Since information about the 
shared triangle edges is already processed and stored, the 
appropriate T-mesh sequence can correspondingly be 
defined by traversing voxel faces and considering rows of 
triangles as a group. A Silicon Graphics VGX or VTX 
hardware can be used to draw surfaces at the rate of more 
than 1.6xl0 6 T-mesh triangles per second. 

In order to further reduce the computational requirements 
of the SpiderWeb method it is possible to simplify the 
computation of the edge hit locations. For example, if an 
edge has a hit, instead of calculating its location exactly, in 
one embodiment of the present invention the hit location is 
taken to a standard position (a preferred choice is midway 
between the edge vertices). The independent angles between 
surface triangles arc thus reduced to a relatively small 
number of cases, which simplifies the display at the expense 
of a somewhat "blocky" appearance of the surface. FIGS. 15 
and 16 illustrate the difference between a fully processed, 
and a blocky approximation surfaces. 

Referring again to FIG. 2, the output structure of the 
SpiderWeb processor 100 can be described as being in a 
vector form, as a sequence of triangle vertex points. In one 
embodiment of the present invention the data supplied to 
display processor 70 comprises a list of such vectors in three 
dimensions. The list also includes the shared normal vector, 
associated with each triangle vertex. Such normal vectors 
are computed at steps 260 and 370 in FIGS. 9 and 12, 
respectively. In a preferred embodiment of the present 
invention, a shaded, rasterized image is produced by the 
display processor 70 and supplied to display 90, in which 
each triangular tile is shaded on a scale, dependent on the 
relation between the viewing angle selected from the control 
50 and the normal vector to the approximation tile. 
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Specifically, in order to produce smoothly shaded images, 
the unit normal of all triangles sharing a common vertex is 
averaged together at the vertex to produce a shared, or vertex 
normal. Each triangle normal is computed as the vector cross 
product of two defining edges. FIG. 17 presents an illustra- 5 
tion of the concepts. The graphics processor 70 can then use 
the shared normals, to grade the calculated intensity over the 
triangle to create a smoothly shaded image having a natural 
appearance. Different coloring can also be used in the 
processor to provide better visual separation of the indi- 
vidual surface tiles, and of different objects or structures. 

Since the computation of the vertex normals in accor- 
dance with the present invention uses localized planar 
triangle information, as opposed to other known methods of 
computing surface normals which typically use the volume 
gradient case, when closely positioned isosurfaces need to 15 
be discriminated, the method of the present invention pro- 
vides a better surface separation. FIGS. 18 and 19 show 
examples of medical imaging surfaces, rendered in accor- 
dance with the method of the present invention. 

It is to be understood that the claimed invention is not to 20 
be limited to the description of the preferred embodiment 
but encompasses other modifications and alterations within 
the spirit of the inventive concept which scope is defined by 
the following claims. 

What is claimed is: 25 

1. A method for rendering a three-dimensional (3D) 
surface structure on a display device of a computer graphics 
display system using information from sequenced image 
slices Produced by a non-invasive data gathering device, 
comprising the steps of: 30 

a) generating a 3D field of data representing the values of 
at least one physical property associated with a 3D 
body by a nonintrusive cross-sectional image gathering 
device; 3J 

b) storing the generated 3D field of data at regularly 
spaced grid locations defining voxels within said body, 
each voxel having eight vertex values and six faces; 

c) selectively retrieving the eight voxel vertex values of 
said physical property for a voxel which is intersected 4Q 
by the surface structure; 

d) determining which voxel vertex values are inside and 
which voxel vertex values are outside of the surface 
structure; 

e) computing, for each voxel face having at least one 45 
voxel vertex inside and at least one voxel vertex outside 
of the surface structure, coordinate values of the inter- 
sect between the voxel face and the surface structure; 

f) forming, for each intersect of a voxel face with the 
surface structure, a triangle having a third vertex in the 50 
interior of the voxel to obtain an approximation of the 
surface structure in the interior of the voxel; and 

g) generating image data associated with the formed 
triangle, said data being displayable at a display posi- 
tion corresponding to the position of the voxel on the 55 
surface structure. 

2. The method of claim 1 wherein step d) further com- 
prises the step of comparing each voxel vertex value with a 
predetermined threshold, the value of said threshold corre- 
sponding to the surface structure. 60 

3. The method of claim 2 wherein at step e) the compu- 
tation of the intersect coordinate values is done by a linear 
interpolation of the voxel vertex values. 

4. The method of claim 3 wherein at step f) the physical 
location of the third vertex of each approximating triangle is 65 
identical and is computed as the cenlroid of all interpolated 
intersect coordinate values. 
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5. The method of claim 4 further comprising the steps of: 

h) identifying each adjacent voxel, sharing a face intersect 
with the surface structure; and 

i) recursively applying steps c), d), e), f), g) and h) to each 
identified voxel until all face-adjacent voxels within the 
three-dimensional body having a face intersect with the 
surface structure are processed. 

6. The method of claim 5 wherein in the step of recur- 
sively computing the intersect coordinates, only coordinate 
values that have not previously been computed are pro- 
cessed. 

7. The method of claim 5 further comprising the step of 
displaying a two-dimensional image of the approximated 
surface structure. 

8. The method of claim 7 wherein the surface structure is 
displayed by utilizing only the identified voxels. 

9. A method for rendering of an internal surface structure 
within a three-dimensional (3D) body on a display device of 
a computer graphics display system using information from 
sequenced image slices produced by a non-invasive data 
gathering device, comprising the steps of: 

a) generating a 3D field of data representing the values of 
at least one physical property associated with the 3D 
body by a nonintrusive cross-sectional image gathering 
device; 

b) storing the generated 3D field of data at regularly 
spaced grid locations defining voxels within said body, 
each voxel having associated with it vertex locations 
and voxel faces, each two intersecting voxel faces 
defining a voxel edge, said voxel edge passing through 
two voxel vertices; 

c) selectively accessing two voxel vertices of a voxel edge 
for a voxel which is intersected by the surface structure; 

d) determining which voxel vertex of the selected voxel 
edge is inside and which voxel vertex of the selected 
edge is outside of the surface structure; 

e) computing, for each selected voxel edge having at least 
one voxel vertex inside and at least one voxel vertex 
outside of the surface structure, a coordinate value of a 
first hit location where the surface structure crosses the 
voxel edge; 

f) examining each of the two intersecting voxel faces 
defining the selected voxel edge for a second hit 
location defining an intersect between the surface struc- 
ture and the voxel face, said intersect passing through 
said first and said second hit locations; 

g) forming, for each intersect of a voxel face with the 
surface of interest, a triangle with a vertex in the 
interior of the voxel to obtain an approximation to the 
surface structure in the interior of the voxel; and 

h) generating image data associated with the formed 
triangle, said data being displayable at a display posi- 
tion corresponding to the position of the voxel on the 
surface structure. 

10. The method of claim 9 wherein step d) further 
comprises the step of comparing each voxel vertex value 
with a predetermined threshold, the value of said threshold 
corresponding to the surface structure. 

11. The method of claim 10 wherein at step e) the 
computation of the intersect coordinate values is done by a 
linear interpolation of the voxel vertex values. 

12. The method of claim 11 wherein at step f) the location 
of the third vertex of each approximating triangle is com- 
puted as the centroid of the interpolated intersect coordinate 
values forming a closed contour around the inside voxel 
vertex. 
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13. The method of claim 12 further comprising the steps 

of: 

i) identifying each adjacent voxel, sharing a face intersect 
with the surface structure; and 

j) recursively applying steps c), d), e), f), g), h) and i) to 
each identified voxel until all face-adjacent voxels 
within the three-dimensional body having a face inter- 
sect with the surface structure are processed. 

14. The method of claim 13 wherein in the step of 
recursively computing the intersect coordinates, only coor- 
dinate values that have not previously been computed are 
processed. 

15. The method of claim 13 further comprising the step of 
displaying a two-dimensional image of the approximated 
surface structure. 

16. The method of claim 15 wherein the surface structure 
is displayed by utilizing only the generated image data 
associated with the formed triangles in the identified voxels. 

17. The method of claim 9 wherein the planar arrange- 
ment of the regularly spaced grid-locations is other than a 
rectangular array. 

18. The method of claim 9 wherein said second hit 
location being closer to the inside voxel vertex as compared 
to the outside voxel vertex of the selected edge. 

19. A computer graphics display system for analysis and 
display of three-dimensional (3D) iso-value surface struc- 
tures using information obtained from sequenced image 
slices produced by a non-invasive data gathering device, 
comprising: 

first storage means for storing three-dimensional signal 
patterns representing the value of at least one physical 
property associated with a three dimensional body at 
regularly spaced grid locations within said body, a 
voxel being defined by a set of adjacent locations in the 
three-dimensions spaced grid, each voxel having voxel 
edges and voxel faces associated with it; 

means for defining a plane of pixel positions having 
preselect position density and corresponding to one 
plane of said spaced grid locations; 

means for determining a pixel on a plane defined by said 
means for defining, said pixel being associated with a 
voxel that is intersected by a surface structure of 
interest within said body; 

processing means for computing coordinate values of the 
cross-section between the surface structure of interest 
and the intersected voxel edges and for computing the 
coordinate values of at least one articulating location 
positioned in the interior of the intersected voxel; 

second storage means for storing the computed coordinate 
values along with connectivity information about join- 
ing the computed coordinate values into a set of 
triangles, each triangle having one triangle edge on a 
voxel face and one triangle vertex at an articulating 
location in the interior of the intersected voxel; 
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display processor means receiving data from said second 
storage means and converting the received data about 
the set of triangles to a display format; and 

means for interactively displaying surface structures of 
interest, said means for interactively displaying being 
driven by said display processor means. 

20. The system of claim 19 wherein said processing 
means further comprises first generating means for gener- 
ating the unit normal vectors of each set of triangles asso- 

' ciated with a location having computed coordinate values. 

21. The system of claim 20 wherein said processing 
means further comprises second generating means for gen- 
erating shared normal vectors associated with the voxel 
locations having computed coordinate values, said shared 

' normal vectors being generated on the basis of the triangle 
normal vectors generated by said first generating means and 
said shared normal vectors being stored in said second 
storage means. 

22. The system of claim 21 wherein the display processor 
3 means uses said stored shared normal vectors for providing 

to the display means image data which, when displayed, 
produces different appearance of different surface structure 
locations being displayed, said appearance being dependent 
on said shared normal vectors. 

5 23. The system of claim 19 wherein the processing means 
further comprises means for identifying which of the voxels 
sharing a voxel face with the intersected voxel are inter- 
sected by the surface structure of interest. 
24. The system of claim 23 wherein the processing means 

3 further comprises means for generating, for each voxel 
identified as being intersected, a plurality of index values 
representing information about the voxel face intersects 
shared, said plurality of index values being stored in said 
second storage means. 

5 25. The system of claim 19 wherein the processing means 
further comprises means for retrieving signal patterns asso- 
ciated with each grid location and means for providing the 
three-dimensional location associated with each said grid 
location. 

0 26. The system of claim 19 wherein said processing 
means comprises means for simultaneously comparing all 
vertex values of said intersected voxel to a value of said at 
least one physical property indicative of the surface structure 
of interest, the output of said means for simultaneously 

5 comparing being used to compute the coordinate values of 
the cross-section independently for each voxel face. 

27. The system of claim 19 further comprising variable 
means for adaptively selecting the position density of the 
pixels of said defined pixel plane. 

0 28. The system of claim 19 further comprising structure 
analysis means receiving input from said second storage 
means and capable of producing on output global structure 
parameters. 

+ + + + * 
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